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Abstract 

Carbonate chemistry variability is often poorly characterized in coastal regions and patterns of covariation with other 
biologically important variables such as temperature, oxygen concentration, and salinity are rarely evaluated. This absence 
of information hampers the design and interpretation of ocean acidification experiments that aim to characterize biological 
responses to future PCO2 levels relative to contemporary conditions. Here, we analyzed a large carbonate chemistry data set 
from Puget Sound, a fjord estuary on the U.S. west coast, and included measurements from three seasons (winter, summer, 
and fall). PCO2 exceeded the 2008-201 1 mean atmospheric level (392 (latm) at all depths and seasons sampled except for 
the near-surface waters {< 10 m) in the summer. Further, undersaturated conditions with respect to the biogenic carbonate 
mineral aragonite were widespread {Q,g^<^). We show that PCO2 values were relatively uniform throughout the water 
column and across regions in winter, enriched in subsurface waters in summer, and in the fall some values exceeded 
2500 (latm in near-surface waters. Carbonate chemistry covaried to differing levels with temperature and oxygen 
depending primarily on season and secondarily on region. Salinity, which varied little (27 to 31), was weakly correlated with 
carbonate chemistry. We illustrate potential high-frequency changes in carbonate chemistry, temperature, and oxygen 
conditions experienced simultaneously by organisms in Puget Sound that undergo diel vertical migrations under present- 
day conditions. We used simple calculations to estimate future PCO2 and values experienced by diel vertical migrators 
based on an increase in atmospheric CO2. Given the potential for non-linear interactions between PCO2 and other abiotic 
variables on physiological and ecological processes, our results provide a basis for identifying control conditions in ocean 
acidification experiments for this region, but also highlight the wide range of carbonate chemistry conditions organisms 
may currently experience in this and similar coastal ecosystems. 
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increased pC02 availability (e.g., photosynthesizers), while others 
(e.g., echinoderms, molluscs, some fishes) may be adversely 
impacted by enhanced shell dissolution rates, reductions in 
metabolic scope, or neurological impairment [2,6-8]. A major 
challenge to inferring the ecological implications of these studies, 
however, is an incomplete understanding of present-day patterns 
of carbonate chemistry variability, particularly in coastal habitats 
[9-14]. 

In continental shelf and estuariiie systems, complex physical 
processes and high rates of biological production act to perturb 



Introduction 

Global surface ocean pH levels have decreased by approxi- 
mately 0. 1 units since the Industrial Revolution and will 
potentially decline an additional 0.3 to 0.4 units over the next 
century due to anthropogenic CO2 emissions [1-4]. This change 
in ocean chemistry, termed "ocean acidification" (OA), results in a 
reduction in the saturation state of carbonate minerals such as 
calcite and aragonite, elevates PCO2, and lowers pH [5]. A 
growing body of research indicates that such changes may impact 
a diversity of species, with some organisms benefiting from 
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pC02 levels away from conditions under air-sea gas equilibrium 
[9] and may exhibit considerable variability over interannual, 
seasonal, and daily time scales [10,14-17]. In contrast, researchers 
studying organisms from these systems often use IPCC projections 
of global average surface ocean carbonate chemistry conditions 
that correspond to assumptions of air-sea gas equilibrium to set 
control and experimental carbonate chemistry conditions [11,18]. 
Recognition of the potential disparity between conditions consid- 
ered "controls" in OA experiments and actual carbonate 
chemistry conditions experienced by organisms in the habitats 
they occupy has been highlighted in a suite of recent publications 
[9,10,18-20]. Consequently, there is increasing awareness that 
researchers should attempt to simulate control conditions that 
correspond to carbonate chemistry conditions experienced by 
organisms in natural habitats and interpret experimental outcomes 
accordingly [21,22]. 

However, largely missing from this discussion is the observation 
that carbonate chemistry commonly covaries with other biolog- 
ically relevant variables including temperature and oxygen. 

In some regions, the strength of these associations has permitted 
the development of semi-mechanistic predictive models of 
carbonate chemistry in coastal environments based on field 
measurements of temperature and oxygen [15,23]. Yet for 
researchers focused on laboratory OA experiments, these 
relationships are rarely incorporated into experimental designs. 
This oversight poses important potential drawbacks. Foremost, if 
carbonate chemistry strongly covaries with temperature and 
oxygen, researchers risk running experiments with control water 
characteristics that are atypical of the habitat that a focal organism 
/ life stage occupies. Physiological parameters including aerobic 
capacity and metabolic scope (the amount of energy that can be 
allocated to activities beyond those for basic existence) are strongly 
influenced by temperature, and adversely impacted by reductions 
in ambient oxygen availability or increases in CO2 concentrations 
[24]. Under extreme warm and cool temperatures defining the 
upper and lower boundaries of the thermal window an organism 
survives within, the negative efiects of oxygen scarcity or high 
concentrations of COg are potentially exacerbated [24] . Because 
of the individual and interactive efiects that temperature, oxygen, 
and PCO2 may have on physiology and other metrics related to 
fitness and survival [6,24,25], knowledge of their covariabUity is 
essential for designing OA experiments that adequately charac- 
terize biological performance under contemporary relative to 
future acidified conditions. 

Here, we have analyzed a large carbonate chemistry data set 
from Puget Sound, a large fjord estuary on the U.S. west coast 
[26]. Puget Sound is a highly altered ecosystem that suffers from 
multiple stressors stemming from extensive urbanization in the 
surrounding watershed. Habitat loss, shoreline development, 
pollution, alteration of freshwater flows, and overfishing are 
pervasive problems [27]. Despite these impacts, the Sound 
remains a biologically diverse system that directly supports 
economically important industries including aquaculture, com- 
mercial and recreational fisheries, and tourism [27]. Located on 
the eastern boundary of the NE Pacific, the Sound is supplied with 
CO2- and nutrient-rich waters that are upwelled onto the 
continental shelf [26,28]. Further CO2 additions resulting from 
high rates of respiration act to make Puget Sound a highly 
vulnerable region for acidification [26]. Due to the reduced 
carbonate buffering capacity of the seawater, increases in pC02 
and reductions in pH resulting from ocean uptake of anthropo- 
genic CO2 will occur more rapidly compared to open-ocean 
waters [26]. Similarly, the saturation state of seawater with respect 
to biogenic calcium carbonate minerals (calcite and aragonite) will 



also decline at a faster rate, suggesting that calcifying species may 
be adversely affected sooner in Puget Sound relative to open ocean 
populations [26]. 

We sought to evaluate seasonal and regional variability in 
carbonate chemistry conditions as well as patterns of covariation 
with temperature, oxygen, and salinity. Our primary rationale was 
to provide an evaluation of carbonate chemistry variability in this 
system as part of a larger effort to develop OA experimental 
designs with improved ecological relevance. Our second objective 
was to evaluate the potential range of carbonate chemistry 
conditions experienced by members of the Puget Sound pelagic 
food web and how that range may change with additional CO2. In 
particular, we focused on potenticd simultaneous changes in 
temperature, oxygen, and carbonate chemistry conditions that are 
likely to be experienced by pelagic biota undergoing diel vertical 
migrations across regions and seasons in Puget Sound. In doing so, 
we aim to draw attention to the more general issue of identifying 
temperature, oxygen, and carbonate chemistry combinations that 
reflect plausible control and experimental conditions that are 
seasonally appropriate and correspond to the habitat occupied by 
the focal study organism or assemblage. We analyzed water 
chemistry data collected over three seasons (winter, summer, and 
fall), which afforded the opportunity to assemble seasonal 
snapshots of spatial variation in carbonate chemistry. Further, 
we collected data over two summers (August 2008 and September 
2009) and two fall periods (October 2010 and October 201 1). The 
inclusion of summer and fall sampling permitted an evaluation of 
interannual variation over a particularly dynamic season in Puget 
Sound. During the transition from late summer to early fall, 
intrusions of denser oceanic waters at depth result in the upward 
movement of the most oxygen-depleted waters in southern Hood 
Canal (a major sub-basin of Puget Sound) and southerly winds 
over Hood Canal result in episodic, localized upwelling events that 
bring oxygen-depleted, cooler waters to the surface at the southern 
end of Hood Canal [29]. We show that carbonate chemistry in 
Puget Sound covaries strongly with temperature and oxygen, but 
that the strength and sign of the relationship can change 
seasonally. Our findings offer insight into patterns of carbonate 
chemistry variability in this system, and we highlight their 
implications for the design and interpretation of OA laboratory 
experiments that simulate conditions representative of this and 
other coastal systems. 

Materials and Methods 

We collected water samples in Puget Sound over five cruises 
that spanned winter (February 2008), summer (August 2008, 
September 2009) and fall (October 2010 and 2011; Table 1) as 
part of the Puget Sound Regional Synthesis Model Program 
(PRISM) conducted by the University of Washington in collab- 
oration with the Washington Department of Ecology. The 
summer surveys correspond to conditions prior to the intrusion 
of denser oceanic waters at depth while the fall surveys capture this 
process. The full survey includes transect lines with stations in all 
four major sub-basins of Puget Sound as well as the Strait of Juan 
de Fuca (Fig. 1). For the current study, we focused on stations 
south of the sill at the entrance to Hood Canal (herein referred to 
as 'Hood Canal') and those spanning the northern entrance of 
Admiralty Inl(;t to northern Central Puget Sound (herein 
'Admiralty Inlet'). The bottom depth of the station near the sill 
in Hood Canal is ~100 m while the deepest stations near the 
center of the basin are ~ 1 70 m. Bottom depths of stations 
spanning Admiralty Inlet range from 70 m to 160 m. We focused 
our analysis on these regions because they possess contrasting 
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Table 1. 


Dates and stations sam 


pled In Hood Canal and 


Admiralty Inlet for carbonate chemistry during PRISM cruises 


from 2008 to 2011. 






Cruise 


Sampling dates 


stations sampled 


Feb. 2008 


5-7 Feb. 


7-22, 24, 401 


Aug. 2008 


12 - 13 Aug. 


7-22, 24, 401 


Sep. 2009 


30 Sep. - 1 Oct. 


7-11, 15-19, 21, 22 


Oct. 2010 


15-30 Oct. 


7-22, 24, 401 


Oct. 201 1 


12 - 13 Oct. 


7-22, 24, 401 


doi:l 0.1 371 /journal.pone.0a8961 9.t001 



mixing pattem.s and were sampled consistently across cruises. AH 
stations were located in Washington State waters and no special 
permission was required to collect water samples. Field activities 
did not involve endangered or protected species. 

At all stations, fuU water column profiles of temperature, 
salinity, and oxygen were obtained using a SeaBird (SBE) 9 1 1 plus 
conductivity-temperature-density (CTD) instrument equipped 
with an SBE 43 oxygen sensor (SeaBird Electronics Inc., Bellevue, 
WA). Data were binned at 0.5-m intervals. To correct oxygen 
sensor measurements for potential offsets, discrete water samples 
were analyzed for oxygen by modified Winkler titration [30] . The 
salinity sensor on the CTD was factory-calibrated annually; the 
sensor was checked against discrete values collected on the cruise 
at a subset of stations, analyzed with a GuUdline Autosal by the 
University of Washington, and calibrated, if necessary. Cruise 
CTD data are available through the NANOOS Visualization 
Server (http:// nvs.nanoos.org/ CruisePrism). 

Discrete water samples were collected and analyzed for 
carbonate chemistry following previously pubfished methods 
[12]. Briefly, water samples were obtained with Niskin-type bottles 
and analyzed in the laboratory for dissolved inorganic carbon 
(DIG), total alkalinity (TA), oxygen, and nutrients. DIC was 
analyzed using coulometric titration [31—33]. TA was measured 
by the potentiometric titration method [33]. Certified Reference 
Materials (CRMs) were analyzed with both the DIC and TA 
samples as an independent verification of instrument calibrations 
[33]. The DIC and TA data have a combined precision and 
accuracy of ~ 1 Jtmol kg~ ' and ~ 2 |J.mol kg~ ' , respectively, based 
on replicate sample analysis and the use of CRMs. Carbonate 
chemistry measurements from all surveys are currenfly being 
prepared for submission and archival at the Carbon Dioxide 
Information Analysis Center (http:/ / cdiac.ornl.gov/ oceans/). 

For our analysis, we focused on carbonate chemistry variation 
as conveyed by changes in pCOj and the saturation state of 
aragonite (flar). We focused on pC02 because this is the most 
commonly used treatment variable in OA experiments and the 
parameter directly changed by anthropogenic COj emissions. We 
evaluated variation in Q.,^^ because of its direct relationship to the 
stability of aragonite material in calcifying species. In supersatu- 
rated conditions (flar >!), the formation of aragonite is 
energetically favored, whUe dissolution is favored in undersaturat- 
ed conditions (£2ar <!)■ We calculated pCOj and Q^r using the R 
library 'seacarb' [34] with dissociation constants from Lueker et al. 
[35]. The pressure effect on the solubility was applied as described 
from the equation of Mucci [36], incorporating the adjustments to 
the constants recommended by MiUero [37]. Based on the 
uncertainties in the DIC and TA measurements, uncertainty in 
the calculated aragonite saturation state and pC02 is approxi- 
mately 0.015-0.04 and 5-10 |J.atm, respectively. 



We evaluated covariation of both pCOa and fi^r with 
temperature, oxygen, and salinity within each survey for Hood 
Canal and Admiralty Inlet stations using simple linear regression 
separately for each variable. We evaluated the degree of 
association based on the variance explained by the regression 
models [R^). To comply with assumptions of normality in the 
response variable, both pC02 and O.^^ were logio-transformed 
prior to fitting the regression models. The fitted regressions offer 
information on linear trends between variables but are not suited 
for the prediction of carbonate chemistry variables using 
temperature, oxygen concentration, or salinity measurements 
beyond the region and time period sampled for each cruise. 

Next, we sought to visually compare seasonal and spatial 
patterns in the vertical distribution of pC02 and flar- To facilitate 
comparability, we predicted pC02 and flar across depths at each 
station using continuous depth profiles of temperature, oxygen, 
and salinity based on relationships estimated with the discrete 
carbonate chemistry samples. The predicted pC02 and Q.^ values 
offered a higher vertical resolution of carbonate chemistry 
variability at each station than measurements from the discrete 
water samples alone. In preliminary stages of model development, 
we first considered multiple regression methods to predict pC02 
and flar, but significant coUinearity in the predictor variables was 
observed in some regions and time periods (absolute value of 
Pearson's correlation coefficient greater than 0.90), which 
hampered their application [38]. To overcome this issue, we 
instead utilized partial least squares regression (PLSR) to model 
pC02 and fiar- PLSR is a dimension reduction technique that is 
weU-suited to data sets that exhibit multicoUinearity among 
predictor variables [39] . Essentially, latent factors (or component 
variables) are obtained from linear combinations of the predictor 
variables. The approach is similar to principal components 
regression except that in PLSR, component variables are obtained 
by maximizing the explained variance in the response variables 
rather than within the predictor variables [39]. We selected the 
optimal number of components by calculating the root mean 
square error of prediction (RMSEP) for each component using 
leave-one-out cross validation. We retained only those components 
without which the RMSEP did not show a significant decrease 
[40]. One PLSR model was estimated for pC02 and a second 
model was estimated for Qa^. We used CTD temperature, salinity 
and oxygen concentration as predictor variables in both models. 
However, to predict relationships for aU stations across time 
periods and regions, we also included a categorical variable that 
denoted the region and cruise from which a sample was obtained 
along with first-order interactions between the continuous CTD 
variables and the categorical variable. Doing so permitted the 
relationships between the carbonate chemistry parameters and 
temperature, salinity and oxygen to vary between regions and 
seasons. This was deemed necessary because preliminary exam- 
ination of the data sets suggested differing relationships between 
the explanatory and response variables across seasons and regions. 
As with the simple linear regression models, both pC02 and fiar 
were logio-transformed to improve normality in the residual error 
structure. To evaluate model adequacy, we evaluated the 
explained variance in the response variable (i?^), RMSEP, and 
checked for homogeneity in the residual error structure. We used 
the fitted model to predict pC02 and Qar values across depths 
using CTD and oxygen sensor measurements and visually 
compared patterns of carbonate chemistry variation between 
regions and seasons; we employed PLSR purely as a means to 
estimate high-resolution vertical profiles of pC02 and Qar- 

We note that although previous studies have developed semi- 
mechanistic predictive multiple regression models of carbonate 



PLOS ONE I www.plosone.org 



3 



February 2014 | Volume 9 | Issue 2 | e89619 



Covariation between CO^ Temperature, and Oxygen 




123°W 123.5°W 122°W 



Figure 1. Map of the study area and stations sampled for carbonate cKiemistry in Hood Canal (denoted by solid circle symbol) and 
Admiralty Inlet (open triangle) from Puget Sound, Washington, USA. Sampling tool< place February and August 2008, September 2009, and 
October 2010 and 2011. Station numbering follows the convention of PRISIM cruises (1998-2011, httpV/nvs.nanoos.org/CruisePrism). 
doi:1 0.1 371 /journal.pone.008961 9.g001 



chemistry using these variables in the California Current system 
[15,23], the application of similar methods to Puget Sound was 
hampered by several considerations beyond the multi-coUinearity 
we observed among predictor variables. Namely, a key assumption 
of these models is that the stoichiometry of the C02:02 
relationship in constituent water masses is primarily controlled 
by aerobic respiration acting over the time span since last exposure 
to the atmosphere. Puget Sound receives waters upwelled onto the 
continental shelf that are entrained in subsurface flows through the 
Strait of Juan de Fuca, but that mix to the surface and combine 
with low salinity, seaward-flowing surface waters due to the 
presence of shallow sills at Admiralty Inlet and the entrance to 
Hood Canal [41]. As a result, COjiOg relationships in the Sound 
reflect a combination of respiration in inflowing upwelled coastal 
waters, differential diffusion rates of CO^ and O2 between 
seawater and the atmosphere, and variable mbdng with surface 
waters where C02:02 relationships are strongly influenced by 
biological production and freshwater inputs. We expect that semi- 
mechanistic models may exist with the capacity to estimate at least 
some carbonate chemistry variables in Puget Sound (e.g. TA), but 
our primary aim here was to interpolate pC02 and il^r values 
between depths where discrete samples were collected at stations 
within each survey. The development of models to extrapolate 
carbonate chemistry parameters across time periods and locations 
outside our spatial and temporal sampling domain is beyond the 
scope of the present study but remains an area of active research. 

In Puget Sound as elsewhere, many fish and zooplankton, 
including euphausiids, mysids, copepods, pteropods, and a large 
variety of fish and invertebrate larvae, undergo diel vertical 
migrations (DVM) and may swim vertical distances in excess of 
50 m in less than an hour [42-46]. In stratified waters, diel vertical 
migrators may pass through water masses with drastically differing 



temperature, oxygen, and carbonate chemistry conditions [45]. 
The prevalence of DVM in Puget Sound fish and invertebrates 
and our need to understand the response of these groups to 
acidification motivated our interest in describing and comparing 
potential carbonate chemistry, temperature, and oxygen condi- 
tions experienced by diel vertical migrators across seasons and 
regions. In particular, we evaluated how these variables may 
change simultaneously as organisms move between near-surface 
and subsurface waters and compared contemporary carbonate 
chemistry conditions experienced by diel vertical migrators with 
those assuming an increase in DIG due to future ocean uptake of 
anthropogenic CO2. 

To estimate pC02 and under elevated DIG conditions, we 
used an equilibrium approach based on the thermodynamic 
erjuations that describe the carbonate system [5] . From the outset, 
we recognized the complexity of physical (vertical and horizontal 
mixing, air-sea CO2 fluxes, freshwater inputs) and biological 
processes (photosynthesis, respiration) that potentially regulate 
carbonate chemistry conditions in Puget Sound and therefore did 
not seek to generate precise predictions of pG02 and Q..^^ 
corresponding to specific IPGG scenarios of GO2 emissions. 
Rather, we considered a scenario in which the present-day 
anthropogenic DIG load in Puget Sound (~30 nmol kg~'; [26]) 
was rjuadrupled to 120 |J,mol kg~'. That is, we added 90 |a,mol 
kg ' to observed in situ DIG measurements from the discrete 
water samples, and recalculated the carbonate system to estimate 
future seawater pG02 and Q^r (hereafter noted as pC02* and 
fiar*)- This DIG addition approximates potential anthropogenic 
CO2 loads in mid-latitude surface ocean waters after equilibration 
with an 788 ^atm pC02 atmosphere as predicted for the year 
2100 based on the Intergovernmental Panel on Climate Change 
(IPCC) IS92a 'continually increasing' emissions scenario [1,5]. We 
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Figure 2. Patterns of covariation between PCO2 and temperature (upper), oxygen (middle), and salinity (lower panels) in Hood 
Canal and Admiralty Inlet, Washington. For reference, the average atmosplneric PCO2 level (2008-201 1 mean: 392 |iatm) is noted by a dashed 
line. Trend lines correspond to relationships fitted using logio-transformed PCO2. 
doi:1 0.1 371 /journal.pone.008961 9.g002 



derived our DIG increment value from potential changes to 
inflowing marine waters to Puget Sound because river discharge 
comprises only 5% of the annual water exchange, and the median 
water residence time is ~80 days for Hood Canal and less than 60 
days for the remaining sub-basins of Puget Sound [47]. 
Consequently, inflowing marine waters largely govern tempera- 
ture, salinity, and nutrient loads in Puget Sound [28,41,48] and set 
initial conditions for carbonate chemistry in bottom waters 
entering Hood Canal as well [26]. The method carries the 
implicit assumption that alkalinity, respiration rates, and circula- 
tion patterns remain invariant. In summer, strong upwelling winds 
on the Washington coast supply deep, dense waters to the Sound 



that are thought to have last been in contact with the atmosphere 
approximately 50 years prior and that acquired their anthropo- 
genic CO2 burden at that time [12]. The anthropogenic CO2 
burden in Puget Sound can therefore lag behind conditions in 
open-ocean surface waters by a few decades, and future DIC loads 
will depend partly on interannual changes to the seasonal 
upwelling regime. Recognizing these complexities, our application 
of a uniform increase in DIC across seasons and locations is meant 
to simplify and facilitate general comparisons of possible future 
carbonate chemistry variability. As with pC02 and O..^^, we used 
CTD temperature, salinity, and oxygen values corresponding to 
PCO2* and Qar* from the discrete water samples to predict pC02* 
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Figure 3. Patterns of covariation between aragonite saturation state (Har) and temperature (upper), oxygen (middle), and salinity 
(lower panels) in Hood Canal and Admiralty Inlet, Washington. Values below 1 (indicated by dashed line) correspond to corrosive conditions, 
those above are supersaturated with respect to aragonite. Trend lines correspond to relationships fitted using logio-transformed Q.^,. 
doi:10.1371/journal.pone.0089619.g003 



and flar* across all depths using PLSR; model selection proceeded 
in a manner identical to that used for pCOg and Q„. 

Based on typical DVM patterns observed in Puget Sound 
zooplankton and fish [42,45], we assumed that zooplankton 
descended from 3 to 50 m at a rate of 1 m min starting at 0600 
hr, and ascended from 50 to 3 m beginning at 1800 hr at the same 
rate. Given this depth pattern over a 24-hour period, we extracted 
the corresponding CTD data and pCOj, pC02*, flar, and flar* 
values predicted across depths from each station in the survey to 
generate time traces. We compared time traces across seasons and 



between Admiralty Inlet and Hood Canal stations to evaluate 
potential high-frequency variation in carbonate chemistry result- 
ing from diel movements, explore how carbonate chemistry 
changes simultaneously with other biologically important vari- 
ables, and examine how carbonate chemistry may differ under 
C02-enriched conditions as experienced by diel vertical migrators. 
To aid comparability, we used the same DVM across seasons and 
locations, but acknowledge that seasonal, regional, and interspe- 
cific variations in DVM may in fact be significant. 
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Figure 4. (a) Predicted versus calculated PCO2 and PCO2* and 
(b) Sl„ and Sl^r* based on partial least squares regression 
models for Puget Sound water samples obtained from 2008 to 

201 1 . pCOj* and Q^,* values correspond to estimates based on an the 
addition of 90 |j.mol kg"' DIC across all samples. 
doi:1 0.1 371/journal.pone.008961 9.g004 

Results 

In Hood Canal, pC02 estimates obtained from discrete water 
samples ranged from 230 to 2610 |iatm, and the majority of water 
samples exceeded atmospheric pC02 levels (mean atmospheric 
pCOa from 2008 to 2011: 389 natm; Fig. 2). Regression intercept 
and slope coefficients and the explained variance (based on the 
logio-transformed response variables) associated with the trend 
lines depicted in Fig. 2 are presented in Table SI. The Hood 
Canal August relationship between pCOg and temperature was 
negative and spanned the largest range of temperatures (8.5-1 7°C) 
relative to the other months surveyed (Fig. 2). In fall, the range of 
temperatures observed in both years was narrower relative to 
summer (~9-12°C; Fig. 2). Temperatures in winter ranged from 8 
to 1 1°C in Hood Canal but showed the opposite relationship with 
pCOj as observed in the summer and fall (Fig. 2). Similar 
temperature-pC02 relationships were evident for Admiralty Inlet, 
but for all cruises the range of temperatures and pC02 levels were 
narrower. In general, variation in pC02 was more tighdy 
associated with temperature in Hood Canal relative to Admiralty 
Inlet based on the higher values observed in Hood Canal across 
aU seasons (Table SI). 

Negative relationships between pC02 and oxygen were evident 
in all seasons and years in both Hood Canal and Admiralty Inlet 
with the exception of Admiralty Inlet in February (Fig. 2, Table 



SI). As with temperature, the larger range of O2 values were found 
in Hood Canal relative to Admiralty Inlet (Fig. 2). In general, the 
ii(:gati\ e relationship between pC02 and O2 was stronger in Hood 
Canal than Admiralty Inlet based on values (Table SI). 
Salinities, in contrast, spanned a relatively narrow range of values 
(~27 to 33; Fig. 2). In both regions, relationships between pC02 
and salinity were positive in August. In the remaining seasons, 
pCOj and salinity relationships were weaker. Overall, pCOa was 
less tightly associated with salinity relative to temperature and 
oxygen across all seasons and both regions (Fig. 2; Table SI). 

As expected, patterns of Q„ were nearly inverse to those 
observed for pC02 (Fig. 3; Table SI). With the exception of 
August, when in situ Q„ values ranged from 0.2 to 2.8, Q^r values 
were ~1 or lower in both Hood Canal and Admiralty Inlet in aU 
remaining months (Fig. 3). Similar to pC02, the range of 
values was highest in Hood Canal and lowest in Admiralty Inlet 
(Fig. 3). Patterns of association between Q..^^ and temperature, 
oxygen, and salinity based on i?^ values were similar to those of 
PCO2 between regions and across seasons (Table SI). 

Partial least square regression models of pC02 and Sl^ 
explained greater than 95 and 94:% of variance in each survey, 
respectively (Fig. 4). RMSEPs (based on the original response 
variable scale) ranged between 4 and 50 |iatm and 0.02 and 0.05, 
for PCO2 and Q^, respectively. The model fits suggested a 
sufficient level of predictive accuracy for estimating values across 
depths based on CTD depth profiles of oxygen, temperature, and 
salinity. We predicted pC02 and Q^r across depths and examined 
temperature and oxygen depth profiles as well, but did not 
evaluate salinity due to the narrow range of values observed in the 
survey. 

Overall, seawater pC02 values were above mean atmospheric 
levels (392 |Xatm) at all depths and across surveys with the 
exception of shallow waters (< 15 m depth) at several stations 
sampled in August 2008 (Fig. 5). Waters were also undersaturated 
with respect to aragonite (n„j<l) at all depths sampled in 
February 2008 and October 2010 and 2011 (Fig. 5). Q^r in near- 
surface waters (< 15 m) at some stations in Hood Canal in August 
2008 and September 2009, however, were supersaturated (fiar 
>1), while several stations at Admiralty Inlet sampled in August 
2008 were supersaturated throughout the water column (Fig. 5). 

In general, large differences in vertical gradients in carbonate 
chemistry, temperature, and oxygen conditions were apparent 
across seasons and regions. In February 2008, CTD depth profiles 
of temperature and oxygen and predicted values of pC02 and fl^r 
indicated waters in Admiralty Inlet were well-mixed, and showed 
littie vertical or inter-station variability (Fig. 5). In Hood Canal, 
several stations were also weU-mixed, though vertical gradients 
were evident at several stations: deep waters were warmer, oxygen 
depleted, rich in pC02, and exhibited lower Sl^^ values relative to 
surface waters (Fig. 5). In August 2008 and Scpt(^mhcr 2009, the 
water column was strongly stratified at several stations in Hood 
Canal, with warm, oxygen-rich, low pC02, supersaturated waters 
occurring in near-surface waters. At depth, the opposite conditions 
prevailed. Variation in water conditions was comparatively lower 
at Admiralty Inlet both across stations and vertically (Fig. 5). 

Stratification was again low at Admiralty Inlet stations in 
October 2010 and 2011 with near-uniform temperature, oxygen, 
and carbonate chemistry conditions throughout the water column 
and little variability between stations. However, in Hood Canal 
oxygen-poor, C02-rich, and low Q.^ waters occurred in the top 
50 m of the water column, while the opposite conditions were 
observed at depth (Fig. 5). Hypoxic conditions (oxygen concen- 
trations less than ~60 |tmol kg ') were observed in waters at 
depths shallower than 20 m at stations in southern Hood Canal in 



PLOS ONE I www.plosone.org 



7 



February 2014 | Volume 9 | Issue 2 | e89619 



Covariation between CO^ Temperature, and Oxygen 




8 10 12 14 16 100 200 300 500 1500 2500 0.5 1.0 1.5 2.0 2.5 
Temperature (°C) O2 (timol kg i) pCOjInatm) Q„ 



Hood Canal ^^^Admiralty Inlet 

Figure 5. Temperature and oxygen CTD depth profiles at stations sampled in Hood Canal and Admiralty Inlet. PCO2 and Oar values are 
model-based estimates obtained using partial least squares regressions. To aid interpretation, vertical dashed lines are overlaid denoting the average 
atmospheric PCO2 (2008 to 2011 mean: 392 jiatm) and the equilibrium aragonite saturation state (Q^r =1)- Profiles are truncated at 100 m to aid 
comparability. 
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October 2011. Across all five surveys, the highest and lowest 
estimated pCOj and flar values, respectively, were observed at 
these same stations and depths: estimates of pC02 exceeded 2400 
and flar values were less 0.4 (Fig. 5). 

Time traces of temperature, oxygen, pC02, and fl^r concen- 
trations corresponding to conditions experienced by hypothetical 
diel vertical migrators are depicted in Fig. 6. As suggested by the 
depth profiles, diel vertical migrators occupying the sampled 
stations likely experience simultaneous and abrupt changes in 
temperature, oxygen, and carbonate chemistry conditions in 
strongly stratified locations, and these were most prevalent in 
Hood Canal and Admiralty Inlet in August 2008 and September 
2009 (Fig. 6). Time traces in August 2008 in Hood Canal, for 
instance, suggest that mean temperature, oxygen concentration, 
n^r, and pCOg experienced by diel vertical migrators in near- 
surface waters (3 m) at night and deep waters (50 m) during the 
day differed by 5.1°C, 214 |xmol kg"', 1.68, and 1100 ^atm, 
respectively (Fig. 6). Importantiy, the ranges across stations in 
temperature, oxygen, and carbonate chemistry conditions in both 
surface and deep waters were also notable during this time period 
in Hood Canal. The larger range in carbonate chemistry variables 
occurred in deep relative to near-surface waters: minimum and 



maximum pCO^ values in near-surface waters were estimated to 
be 250 and 480 |xatm and in deep waters values were 980 and 
2030 |J.atm, respectively. Qar ranged from 1.52 to 2.68 in near- 
surface waters and 0.37 to 0.77 in deep waters (Fig. 6). 

In all seasons and regions, pC02* values were not only higher 
relative to PCO2, but the difference between near-surface and 
deep-water values increased (Fig. 6). In Hood Canal August 2008, 
the mean difference between near-surface and deep-water pC02 
widened from 1100 to 2090 [Xatm with the uniform addition of 
Die (Fig. 6). Correspondingly, fl^r decreased across all depths as 
well, but the mean difference between near-surface and deep 
water was lower (1.68 vs. 0.98) relative to fiar, as the decline in the 
carbonate ion concentration slows as it approaches zero (Fig. 6). 

In contrast to August 2008 and September 2009, the mean 
differences in temperature, oxygen, and carbonate chemistry 
variables between near-surface and deep waters were lower in 
October 2010 and 2011 (Fig. 6). However, the upward displace- 
ment of the waters with the lowest oxygen concentrations toward 
the surface in Hood Canal resulted in time traces in which diel 
vertical migrators encountered higher pC02 and lower flar during 
the night rather than during the day (Fig. 6). Further, the range in 
PCO2 and flar conditions in Hood Canal was larger among 
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Figure 6. Time traces of temperature, oxygen, PCO2 and fl^r conditions potentially experienced by hypothetical diel vertical 
migrators at stations sampled in Puget Sound, Washington, U.S.A. Organisms ascend from 50 to 3 m at a rate of 1 m min"' beginning at 
1800 hr and descend back to 50 m at 0600 hr. Colored bands correspond to maximum and minimum PCO2 and Qar values and PCO2* and Q^* as 
estimated assuming a 90 |imol kg^' increase in DIC in Hood Canal and Admiralty Inlet stations. The solid lines correspond to mean values. For 
comparison, atmospheric PCO2 (2008-201 1 average: 392 |j.atm) is indicated by a dashed line. Qa, values below 1 (denoted by dashed line) indicate 
undersaturated conditions with respect to aragonite. 
doi:10.1371/journal.pone.0089619.g006 



stations in near-surface relative to deep waters, and these ranges 
increased and decreased, respectively, with an increase in DIC 
(Fig. 6). 

Overall, the difierence in near-surface and deep temperatures, 
oxygen concentrations, and carbonate chemistry conditions 
experienced by diel vertical migrators in February 2008 were 
lower than in the remaining surveys (Fig. 6). Similar to August 
2008 and September 2009, oxygen concentrations and Q,^^ were 
lower and pC02 higher for diel vertical migrators at depth during 
the day in Hood Canal, but unlike the remaining months surveyed 
the mean temperature was warmer by 1.7°C relative to near- 
surface waters experienced during the night. Further, the potential 
range of water conditions experienced by diel vertical migrators at 
a given depth across stations was narrower than observed in the 
remaining surveys (Fig. 6). 

Discussion 

Our study expands on previous research on Puget Sound 
carbonate chemistry that noted waters gready enriched in CO2 
relative to air-sea equilibration conditions in the winter and 
summer [26]. Here, we show that C02-enriched conditions persist 
through the fall. We also show that carbonate chemistry conditions 



covary with temperature and oxygen but that the direction and 
strength of these correlations change with season and between 
nearby regions characterized by differing levels of water mass 
mixing. Our findings suggest that fish and invertebrates that 
exhibit DVM behavior in stratified locations are likely to 
experience large, rapid changes in ambient carbonate chemistry 
conditions over the course of a single day during the summer and 
early fall months. Further, we show through simple calculations 
that a uniform increase in DIC not only elevates pC02 and lowers 
flar, but also widens the absolute difference in these values 
between near-surface and subsurface waters. Together, our 
findings highlight the physicochemical complexity of coastal 
systems and the importance such data hold for the design of 
ecologically relevant OA experiments. 

Our observations in September 2009 and October 2010 and 
2011 add an additional layer of information to our general 
understanding of seasonal variability in carbonate chemistry in 
Hood Canal. Due to sluggish circulation and the availability of 
nutrients in inflowing marine and terrestrial waters, hypoxic 
conditions commonly occur in Hood Canal during the summer 
months [49]. As evidenced by oxygen and pC02 depth profiles 
in August 2008 and September 2009, aerobic respiration 
simultaneously draws down O2 at depth and increases pC02. 
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However, we show that bottom-water intrusions into Puget Sound 
and Hood Canal during October 2010 and 2011 resulted in the 
upward movement of low-oxygen waters as observed previously 
[45] but that these near-surface waters were C02-replete and 
exhibited the most extreme pC02 and Cl^ values observed across 
all surveys. Fall bottom-water intrusions into the Sound occur 
during neap tides and their onset is influenced by storm-related 
fluctuations in the horizontal density gradient spanning the 
double-silled entrance at Admiralty Inlet [50,51]. Although such 
intrusions play a critical role in the renewal of deep water layers in 
Puget Sound, they also result in the surfacing of low oxygen, high 
pC02 waters which, if upweUed suddenly to the surface due to 
favorable southerly winds, can result in extensive fish kiUs [29]. 
Our observations suggest that organisms that reside in shallow 
waters to avoid low-oxygen, C02-rich waters at depth in the 
summer may be unable to avoid exposure at some locations in 
Hood Canal as these waters circulate upward in the fall. 
Combined, our carbonate chemistry data set indicates that 
pC02 values are spatially and seasonally heterogeneous in Puget 
Sound as in other estuarine systems [52] and further underscores 
the need for researchers to consider natural variability in the 
design and interpretation of OA experiments. 

Importandy, we show that carbonate chemistry conditions also 
covary with temperature and oxygen but that the strength and 
slope of the relationships change with season. Failure to consider 
covariation in carbonate chemistry with temperature or oxygen 
may result in experimental designs with 'control' treatments that 
have little correspondence to conditions tliat animals may have 
adapted to in nature and 'acidified' treatments that do not account 
for naturally high pC02 levels in many estuarine and coastal 
waters. For instance, strong negative relationships between pC02 
and temperature in Hood Canal for August 2008 suggest that 
appropriate control pC02 values for experiments on a widely 
distributed animal from this region would depend on the 
temperatures considered in the experiment. Temperatures of 10 
and 15°C occur within the top 50 m of the water column and have 
corresponding control pC02 values of 1270 and 360 ^atm, 
respectively, based on the linear regression between pC02 and 
temperature. Assuming a uniform increase in 90 |J,mol kg ' in 
Die, the respective acidified treatments are 2150 and 760 |iatm, 
again based on the regression between pC02* and temperature. In 
contrast, experimental controls based on winter Hood Canal 
conditions that use temperatures of 8 and 10°C would have pC02 
values corresponding to 850 and 1360 and pC02* values of 1380 
and 1930 |tatm, respectively. Although there is increasing 
recognition that researchers should consider control treatments 
that reflect carbonate chemistry conditions experienced by the 
organism under study [1 1,18,21], little attention has been given to 
understanding patterns of carbonate chemistry covariability with 
temperature and oxygen. Given the combined elfect environmen- 
tal temperature, oxygen, and pC02 have on metabolic scope [53] 
and evidence indicating non-additive, interactive effects of these 
variables in ecological studies [6,54], we recommend care be taken 
in selecting appropriate values that correspond to those likely to be 
experienced by organisms in the field. Doing so may facilitate 
inferences from experimental results to population-level effects. 

High-resolution carbonate chemistry time series have revealed 
considerable diel variability in coral reef [55,56], kelp forest 
[10,57], coastal upwelling [58], and shallow estuarine systems [59] 
that is primarily driven by oscillations in photosynthesis and 
respiration rates. For animals that undergo DVM in vertically 
heterogeneous regions, however, exposure to water masses with 
differing carbonate chemistry, temperature, and oxygen charac- 
teristics is also possible but has received only limited consideration 



in the context of OA [60]. In Puget Sound, DVM behavior is 
common in a diversity of fish and invertebrate species, including 
zooplankton occurring in stratified regions of Hood Canal [42,45]. 
Developing experiments with temporally dynamic carbonate 
chemistry, temperature, and oxygen conditions that mimic the 
experience of organisms in the wild may offer more accurate 
baselines upon which responses to acidified treatments can be 
compared. Our hypothetical time traces offer a basis for 
developing such experimental designs and offer important insights. 
Foremost, we show that during the summer in particular, 
temperature, oxygen, and carbonate chemistry conditions simul- 
taneously covary with depth, a consideration that has previously 
been overlooked in OA experiments focused on diel migrating 
animals from stratified waters [60]. Organismal performance (e.g., 
growth) may potentially improve in dynamic pC02 treatments 
that emulate natural diel oscillations in carbonate chemistry 
relative to static pC02 treatments [61], but we are unaware of 
studies that have simultaneously varied pC02, temperature, and 
oxygen to replicate the fuU suite of environmental changes 
associated with DVM. The technical challenges associated with 
implementing dynamic, multi-stressor experiments are consider- 
able but OA experimental facilities are becoming increasingly 
sophisticated and such systems are currentiy under development at 
several academic and government laboratories. 

Our recalculation of the carbonate system to generate estimates 
of potential future pC02 and fl^r values experienced by animals 
undergoing DVM showed that pC02 and fl^r increased and 
decreased, respectively, at all dc-pths as expected with an increased 
CO2 load. However, the potential difference in pC02 values 
between near-surface and deep waters increased relative to 
contemporary observations, while the potential surface-to-deep 
difference for Cl^^ decreased. This pattern results from nonlinearity 
in the carbonate system at near-surface ocean pH values: as DIG is 
added in the form of CO2 (either from anthropogenic CO2 
emissions or respiration), the carbonate ion concentration of 
seawater declines, while the CO2 concentration and pC02 forms 
of inorganic carbon increase exponentially [62,63]. In waters 
already burdened with high concentrations of CO2 due to 
respiration, an increase in DIC results in an increase (decrease) 
in PCO2 (fiar) that will be larger (smaller) relative to less burdened 
waters [62,64]. Our estimates of potential change in carbonate 
chemistry variability as experienced by DVM fauna are qualita- 
tively similar to simulations in other studies that have shown an 
amplification of natural variability in carbonate chemistry 
associated with increased DIC [65]. These findings suggest more 
broadly that changes in the mean as well as the variance of pC02 
will accompany increased DIC loadings. The influence of changes 
to pC02 variance on organisms is poorly understood and warrants 
further attention. 

In addition, we show that although spatial and vertical variation 
in carbonate chemistry is considerable in Puget Sound, additional 
DIC may place pC02 and fi^r at most stations and depths outside 
the range of values observed presenffy across depths and stations. 
This possibility was most evident in winter where there was no 
overlap between present-day and future pC02 and Qar as 
experienced by fauna undergoing DVM. During winter, the water 
column is less stratified and conditions are more uniform across 
stations and depths [41]. Consequently, lowpC02, supersaturated 
refugia are absent and carbonate chemistry is relatively homog- 
enous. In contrast, some near-surface waters in August remained 
supersaturated (Qar) and the pC02 at levels low, despite the DIC 
addition. These same near-surface waters initially possessed low 
DIC concentrations due to biological drawdown through photo- 
synthesis and thus higher buffering capacity. 
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The general seasonal patterns observed suggest that species with 
COg-sensitive life history stages that occur in winter are even less 
likely to find depths or regions with favorably low pCOj conditions 
in the future. For species with C02-sensitive life history stages that 
occur in summer, the potential for carbonate chemistry refugia 
persisting in the future may be better because high rates of primary 
production can locally drive pC02 levels below air-sea equilibrium 
in the surface layer, though the benefits of such refugia depend on 
the suitability of other biological variables and their persistence in 
time and space. Our analysis also suggests that COj-sensitive 
species may be more fully excluded from regions such as 
Admiralty Inlet that are well-mixed and therefore lack low 
pCOj refugia at most depths and stations relative to more 
heterogeneous regions like Hood Canal. 

We limited our ocean acidification scenario to a uniform change 
in DIG. Our estimates therefore do not account for processes that 
may further alter carbonate chemistry and the vertical and spatial 
distribution of CO2, which may include climate variables 
(upweUing winds, circulation patterns, river discharge), changes 
in watershed use, eutrophication [9], and atmospheric deposition 
of nitrogen and sulfur from shipping vessels [66] or urban 
pollution [67]. Processes that reduce salinity (and thus TA) or that 
increase respiration (and DIG) increase pG02 [64], while warming 
directly affects carbonate system thermodynamics, resulting in 
solubility-related increases in PGO2, but more importandy 
influences ecosystem metabolism, which strongly affects carbonate 
chemistry [68]. Biogeochemical models of coastal and estuarine 
systems are increasingly utilized to identify the potential processes 
that dampen or exacerbate acidification [69,70] and may help 
clarify the contribution of climate variability and different human 
activities to future changes in Puget Sound's carbonate chemistry. 
We also note that our data set does not permit an evaluation of 
diurnal carbonate chemistry patterns, which can be pronounced in 
coastal surface waters, and is limited to stations in open-water 
habitats. We still lack information on shallower habitats such as 
mudflats, eelgrass meadows, pocket estuaries, and oyster beds, 
which are important components of the Puget Sound ecosystem. 
We recommend targeted sampling in these key habitats over 
seasonal cycles and the development of long-term monitoring 
programs in order to better identify baseline carbonate chemistry 
variability and to facilitate the attribution of the contribution of 
anthropogenic GO2. 

Puget Sound receives waters upwelled along the Pacific Goast 
that are already burdened with anthropogenic CO2 [26], and 
concentrations will continue to increase over the proceeding 
decades regardless of possible changes to atmospheric pG02. This 
is because upwelled waters were last in contact with the 
atmosphere ~50 years prior to surfacing and therefore carry 
GO2 burdens that reflect air-sea equilibration at that time [12]. 
Gonsequently, waters already in transit to upweUing centers on the 
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